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Determining the location of the GeV emission in powerful blazars 
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An issue currently under debate in the literature is how far from the black hole is the Fermi-observed GeV 
emission of powerful blazars emitted. Here we present a diagnostic tool for testing whether the GeV emission 
site is located within the sub-pc broad emission line (BLR) region or further out in the pc scale molecular torus 
(MT) environment. Within the BLR the scattering takes place at the onset of the Klein-Nishina regime, causing 
the electron cooling time to become almost energy independent and as a result, the variation of high-energy 
emission is expected to be achromatic. Contrarily, if the emission site is located outside the BLR, the expected 
GeV variability is energy-dependent and with amplitude increasing with energy. We demonstrate this using 
time-dependent numerical simulations of blazar variability and discuss the applicability of our method. 



1. Introduction 



Blazars are by far the most common objects de- 
tected in the gamma-ray sky pQ. Fermi has detected 
blazar variability as short as a few hours [e.g. [3]. Dur- 
ing these flares, the GeV luminosity has been known 
to increase by a factor of up to several compared to 
its pre-flare luminosity. Because blazars cannot be re- 
solved at these energies (or at any other energy, with 
the possible exception of VLB A observations), it is 
impossible to determine the location of these flares by 
direct detection. To address this issue, we propose 
a diagnostic test that utilizes Fermi variability data 
of short flares to determine the location of the GeV 
emission in high-power blazars, namely flat spectrum 
radio quasars (FSRQ). 
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T- 2 AttR 2 c 
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In this work, all primed quantities refer to the co- 
moving frame of the blob and unprimed quantities 
refer to the galaxy frame. If we assume a nominal 
FSRQ accretion disk luminosity of Ldisk ~ 10 45 erg 
s" 1 and that a fraction £ = 0.1 of this radiation is 
reprocessed by both the BLR and the MT, a typical 
luminosity of the external radiation field is L ext ~ 
10 44 erg s" 1 [12]. 

If the emission site is located within the BLR (at 
R ~ 10 17 cm),the photon field can be considered 
isotropic in the galaxy frame and using Eq. [T] its co- 
moving energy density is 



2. Sources of Seed Photons 
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Relativistic effects determine which photon field is 
dominant at varying distances from the central black 
hole, and as a result the location of the GeV flaring 
site determines the dominant source of seed photons. 
The co-moving (jet frame) energy density of a radi- 
ation field V scales as differing factors of T depend- 
ing on the direction from which the photons enter the 
emitting region [SJ [TU] . If the photon field is isotropic 



U' « T 2 U 
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For photons entering the emitting region from be- 
hind the relativistically moving blob 



where T 10 = 10 _1 r, L B lr,44 = 10 _ L BLR , and 
Rblrm = 10~ 17 Rblr- Similarly, the MT photon 
field is isotropic inside the BLR and its co-moving 
seed photon energy density is 
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where L 



MT,44 



10" 44 L 



MT 



and R 



MT, IS, 



10 1s Rmt- Clearly, inside the Rblr the co-moving 
BLR photon field energy density dominates over that 
of the MT by a factor of ~ 100. 

If the emission site is located at R > 10 18 cm (within 
the MT) then the BLR UV photons enter the emitting 
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region practically from behind, and as per Eq. [2] 

U' BLR « 2.6 x 10~ 4 L B lrm r ~mt is t w er § cm_3 - 

(5) 

The IR photons from the MT retain the same co- 
moving energy density previously given by Eq. [4] In 
this case, therefore, it is the MT that dominates the 
co-moving photon energy density. 

These external photon field co-moving energy den- 
sities need to be compared to the synchrotron photon 
field energy density. If Ruob is the size of the emit- 
ting blob, the co-moving synchrotron photon energy 
density is U' synch L synch / '(47rci? 2 , ob r 4 ). The most 
plausible assumption for the size of the emitting re- 
gion, however, is to set an upper limit to it by its 
variability timescale: Rbiob — ct var 8. We then obtain 
a lower limit for the synchrotron energy density 

U'synch ~ 2 -3 L synchA6 t^ar.lh T W er S Cm ~ 3 > ( 6 ) 

where L synchAe = L s x 1CT 46 , t var ^ h = 3600 x t var 
is the observed variation time in hours, and Tio is 
the same as previously defined. Note that a 6 hour 
variability scale has been observed by Fermi [3J. 

Comparison of the external energy density U ex t (us- 
ing Eqs. [l] or [2] depending on if the emission is 
isotropic) and synchrotron seed photon energy den- 
sity U sync h (Eq. [6| shows that an isotropic external 
photon field dominates, i.e U ex t > U sync h, if 

( J P2 \ V* 

T , „ / ^ synch, 46 ^ext, 18 \ 
T > 17.5 —2 (7) 

\ l var,h^ext,44 ) 

Using the BLR and MT descriptions adopted above 
and t var = 6 hours, results in a T > 6.2 inside the 
BLR and T > 11.2 inside the MT, but outside the 
BLR. VLBI studies of superluminal speeds in powerful 
FSRQs [e.g. see figure 24 of IJorstad et al. 2005] show 
that for most FSRQs 10 < T < 20. Therefore, this 
condition is clearly satisfied inside the BLR, and it is 
expected to be satisfied in the majority of cases in the 
MT. 



3. Cooling in the BLR vs Cooling in the 
MT 

The critical difference between the BLR and the 
MT is the energy of the seed photons: photons orig- 
inating from the BLR are UV photons (cq « 10~ 5 J 
while photons originating from the MT are IR pho- 
tons (cq £3 W~ 7 ), where eo = hv/mc 2 . This difference 
by a factor of ~ 10 2 in typical photon energy is critical 
in that it affects the energy regime in which electron 
cooling takes place, and thus the energy dependence 
of the electron cooling time. 



In powerful FSRQs the IC emission in high states 
can dominate over the synchrotron emission by a fac- 
tor of up to ~ 100 [2]. In cases of high Compton 
dominance, the primary electron cooling mechanism 
is IC scattering. For electrons cooling in the Thom- 
son regime (760 ^ 1, where both the electron Lorentz 
factor 7 and the seed photon energy e are measured 
in the same frame) the cooling rate 7 cx 7 2 . For elec- 
trons with 760 > 1 cooling takes place in the Klein- 
Nishina (KN) regime with 7 oc In 7 [5]. The electron 
cooling time given by 7/7 is then cx 7 -1 in the Thom- 
son regime and cx 7/ In 7 in the deep KN regime. In 
the transition from Thomson to KN regimes (around 
760 ~ 1) the cooling time flattens and becomes essen- 
tially energy independent (Fig [I]). 

The effects of the transition between Thomson and 
KN regimes on the electron energy distribution (EED) 
and the resultant spectrum of the synchrotron and IC 
emission have been studied before [e.g. QTJ [HI 
OH OH HE1 CEH1 OH]. In short, because 7 cx 7 2 in the 
Thomson regime and 7 cx In 7 in the KN regime, the 
cooling time t' coo1 =7/7 scales as 7 -1 in the Thomson 
regime and as 7/ In 7 in the KN regime. 

Because the cooling time is approximately energy 
independent around 760 ~ 1, this energy-independent 
cooling time will be manifested at energies lower by 
a factor of ~ 100 for cooling taking place in the BLR 
compared to cooling taking place in the MT, since the 
BLR seed photons have an energy higher than that of 
the MT by a factor of ~ 100. This can be seen in Fig. 
[l] where we plot the electron cooling time for a source 
with a ratio of external photon field energy density U' Q 
in the co-moving frame to co-moving magnetic field 
energy density Ub, U' /Ub = 100. This corresponds 
to a factor of ~ 100 Compton dominance (ratio of 
inverse Compton to Synchrotron luminosity), similar 
to what is observed in the most Compton dominated 
sources. 

The transition from Thomson cooling to KN cool- 
ing, and from KN cooling to synchrotron cooling can 
also be seen in Fig. [T] If the electrons are cooling 
on photons from the BLR, cooling takes place in the 
Klein-Nishina regime and the cooling time scale is ap- 
proximately energy independent around 760 ~ 1 (Fig. 
[l] top panel) . If the electron population cools on pho- 
tons from the molecular torus, cooling takes place in 
the Thomson regime (Fig. [I] bottom panel). The 
cooling time is heavily energy dependent, and any 
variations should consequently exhibit heavy energy- 
dependence. 



4. The Diagnostic Test 

The energy dependence of the cooling time results 
in an energy dependence (or lack thereof) of varia- 
tions: if the blazar emission site is located within the 
BLR, variations should be achromatic. The energy 
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Figure 1: Cooling time in the blob frame as a function of 
co-moving electron energy 7. Top panel: Blazar emission 
site located in the BLR. Bottom panel: Blazar emission 
site located in the MT. The dotted lines represent the 
various cooling mechanisms (blue = Thomson cooling, 
red = KN cooling, green = synchrotron cooling), the solid 
black line is the total cooling time. The dotted black line 
in the top panel indicates where 760, blb ~ 1. Plots were 
calculated for the following values: seed photon energies: 



to,BLR = 3 x 10 and e 



6 x 10 , energy 



densities Ublr = 2.65 x 10~ 2 erg cm" 3 , Umt = 3 x 10~ 4 
erg cm -3 , a Ub that corresponded to a Compton 
dominance of ~ 100, and Ybuik = 10. Because BLR 
photons have energies higher than MT photons by a 
factor of ~ 100, transition from the Thomson to the KN 
regime occurs at values of 7 ~ 100 times lower in the 
BLR (top panel) than in the MT (bottom panel). 



dependence of the variations can be used as a diag- 
nostic test to determine if the GeV emission site is 
located within the BLR. By comparing Fermi light 
curves of flares at different energies, we propose that 
the energy dependence of the light curve can be used 
as a diagnostic test to rule out whether the GeV flare 
originates in the BLR or MT. 



4.1. Numerical Simulation Results 

To demonstrate the effect of the energy indepen- 
dence or dependence of the electron cooling time on 
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Figure 2: Light curve at various energies (BLR seed 
photons): e = 100 MeV, e = 1 GeV , e = 10 GeV. 



the variability of a flare, we utilized a one-zone nu- 
merical model to simulate a flare. We initialized the 
code with values appropriate for a high power blazar 
with a Compton dominance of ~ 100. For this partic- 
ular simulation we assumed a source size R = 10 16 cm, 
bulk Lorentz factor V — 10, co-moving injected elec- 
tron luminosity L — 2 x 10 44 erg s _1 , maximum 
electron Lorentz factor ~f max — 10 5 , and electron in- 
dex p = 2.5. For the case of a flaring region located 
within the BLR we assumed an initial photon energy 
£o = 3 x 10~ 5 and an energy density (in the galaxy 
frame) Ublr = 2.6 x 10 -2 erg cm -3 . For the case of 
a flaring region located outside the BLR we assumed 
an initial photon energy eo — 6 x 10~ 7 and an energy 
density Umt = 3 x 10~ 4 erg cm -3 . For each case 
the magnetic field B was fixed to assume a Compton 
dominance Uec/Ub = 100. 

The system behaves as expected, showing a notice- 
able difference in the decay rate as well as the am- 
plitude of the flare depending on if the seed photons 
originated from within or from outside the BLR (see 
Figs. [2] and [3]). This difference in decay rate and 
amplitude can be used as a diagnostic test to differen- 
tiate between flares that take place inside or outside 
the BLR. 



4.2. Light-Crossing Time Effects 

Because the GeV emitting region is not a point 
source, any change in the luminosity of the blob will 
not be seen instantaneously. Instead, the observed de- 
cay time of the light curve is the result of the actual 
decay time and the light-crossing time inherent in the 
blob. To test whether the light-crossing time would 
erase any difference in the light curves (as predicted 
by our diagnostic) for the case of the emitting region 
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Figure 3: Light curve at various energies (MT seed 
photons): e = 100 MeV, e = 1 GeV, e = 10 GeV. 




Figure 5: Upper limit in the MT energy density (U) as a 
function of bulk Lorentz factor T plotted for At 
corresponding to 5, 10, and 15 h. Accepted values of of U 
and r values fall below the curves. 




Figure 4: Light curves with light-crossing time effects 
factored in. Initial flux and light-crossing time are 
normalized to 1. Light curves are plotted for energies 
differing by a factor of 10. 



being located outside the BLR, we modeled the flux 
of a source with light-crossing time t^c that is de- 
creasing in flux. We assumed the flux of the source is 
decaying exponentially, F(t) = F a exp~*/* c , where t c 
is the cooling time at a specific energy. 

Our diagnostic predicts differences in the decay time 
of the light curves for the case where the emitting 
region is located within the MT; for the purposes of 
this demonstration we assume cooling occurs in the 
Thomson regime and as a result t c cx eT 1 ! 2 . We plot 
the resultant light curves for cooling times differing by 
a factor of 1/^/10 (i.e. energies differing by a factor 
of 10). As evident in Fig. |4j even with light-travel 
time effects convolved with exponential decay times, 
the differences in the decay times are still manifested. 



4.3. Cooling time difference constraints 
on U M t and r. 

A flare that occurs in the MT should exhibit en- 
ergy dependent cooling times because cooling occurs 
in the Thomson regime. Practical application of the 
diagnostic hinges on the requirement that the decay 
times in different energy bands should have detectable 
differences. In the Thomon regime, the cooling time is 
given by t c — frac2>m e cAUT 2 aT e f 1 ^' 2 ', where eo is the 
seed photon energy and e is the photon energy after 
IC scattering. The expected difference in cooling time 
at two different energies is then 



At 



3m,,ce, 



1/2 



AUT 2 a T 



c-V2 

~LE 



,-1/2 
-HE 



(8) 



where €le and cre are the photon energies at which 
we calculate the cooling time, and U is the MT energy 
density. Because of Fermi 's three hour time resolution, 
At > 3 h and this sets an upper limit to the product 
UT 2 . The T-U parameter space for which the dif- 
ference in decay times are distinguishable is plotted 
in Fig. [5] for three plausible values of At, with the 
area below each line corresponding to the acceptable 
part of the parameter space. The permitted part of 
this diagram comfortably accommodates bulk Lorentz 
factors between 10 and 20 and MT photon densities 
10 -4 erg cm~ 3 



4.4. SSC in the MT? 

Our discussion following equation Q suggests that 
it is possible that in the MT, for relatively low val- 
ues of F, U synch > U ec and the dominant pho- 
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ton field is from jet synchrotron photons. Even if 
this is the case, cooling still occurs in the Thom- 
son regime and the energy dependent cooling time is 
still manifested. For FSRQs, the peak of the syn- 
chrotron spectral energy distribution is at v s ~ 10 13 
Hz |13j . which in the comoving frame corresponds to 
e s = hv s /{m e c 2 Y) s» lCT 7 /r. These seed photons, 
therefore, have a lower energy than the energy of the 
MT photons in the comoving frame, and will be deeper 
in the Thomson regime compared to the MT photons. 
Because cooling still occurs in the Thomson regime, 
the energy dependent cooling time, therefore, is ex- 
pected to be present, even in the case that the dom- 
inant photon density is that of the synchrotron pho- 
tons. 



5. Conclusions 

We have presented a diagnostic test that utilizes 
blazar variability to determine the location of the GeV 
emitting site in blazars. The energy difference in seed 
photons originating from the BLR versus seed photons 
originating from the MT causes electrons within the 
emitting site to cool in different energy regimes. 

For the case where the GeV emitting site is located 
within the BLR, cooling takes place at the onset of the 
KN regime, and the resultant electron cooling time 
is energy- independent. We have demonstrated that 
the associated light curves exhibits decay times that 
are approximately energy independent. Conversely, 
for the case where the GeV emitting site is located 
outside the BLR, cooling takes place in the Thom- 
son regime and the electron cooling times are heavily 
energy dependent. In this case, the associated light 
curves exhibit energy dependence of their decay times. 

The energy dependence of the decay time of the 
light curves is visible within the Fermi energies; these 
differences can be used as a diagnostic test to deter- 
mine whether the GeV emitting region is located in- 
side or outside the BLR. These effects are observable 
within the time resolution of ~ few hours that Fermi 
has achieved for bright flares and are not erased due 
to light-travel time effects. Light curves from a suffi- 
ciently bright and rapid flare [such as that in 3C 454.3; 
H] should be compared at different energies. If the 
GeV emitting site is located within the BLR, the de- 
cay times will exhibit no energy dependence, whereas 



if the emitting site is located within the MT, the decay 
times will exhibit energy dependence. 
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